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ABSTRACT 


The oxidation of nickel by atomic oxygen at pressures from 1 to 45 Nm“^ 
between 1050 and 1250®K has been investigated. In these ranges, the oxidation, 
was found to follow the parabolic rate . law, viz. , kp = 1 . 14xl0r^ expC“13410/T).g2cm~ 
sec~^ for films of greater than 1 pm thickness and was pressure independent. 

The activation enthalpy for the oxidation reaction was.112 ±l| kJ mole"^ (27 ±3kcal 
mole"^). Of a number of, possible mechanisms and defect structures considered, 
it was shown that the most likely was a saturated surface defect model for 
atomic oxidation, based on reaction activation enthalpies, impurity effects, . 
pressure independence, and magnitudes of the rates. A model judged somewhat., 
less likely was one having doubly ionized cationic defects rate controlling in 
both atomic and molecular oxygen. From comparisons of the appropriate proces- 
ses, the following enthalpy values were derived: AH* . (Ni diffusion in 

NiO) = 110 ±30 kJ mole"^ and AHf (doubly ionized cation vacancies in NiO from 
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atomie-sbxygen ;= ,,-.9 >25 kJ mole"^ 
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SYMBOLS 


C, C, C" 
D 

(d) 

AG° 

AH° 

AH* 

K 

k, k' 

(n) 

P 

•R 

AOO. 

AS 

(s) 

(s d) 

T 

t 

w . 

X 


y 

[ ] 
□ 


collective constants 

diffusion coefficient 

doubly ionized vacancy 

standard free energy of formation 

standard enthalpy change for reaction 

enthalpy of activation 

equilibrium constant 

rate constant 

non- ionized vacancy 

pressure 

gas constant 

standard entropy change for reaction 
singly ionized vacancy 

singly to doubly ionized vacancy reaction 

absolute temperature 

time 

weight ■ . 

thickness . 

defect mole fraction 
concentration 
vacancy 


Subscripts 

a, c refer to reactions (a) and (c) , respectively 

vacancy diffusion 


diff 


f 


formation 


Ni,0,02 refers to process or constant involving particular element 
o initial state 

p parabolic 

a conductance 
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RATES AND MECHANISMS OF THE ATOMIC OXYGEN REACTION 
WITH NICKEL AT ELEVATED TEMPERATURES 
Jerry D. Christian* and William P. Gilbreath 
Ames Research Center, NASA, Moffett Field, Calif., 94035 


The oxidation rates of nickel and the defect structure and transport 
properties of nickel oxide in molecular oxygen have been the subjects of num- 
erous investigations. It is generally agreed that NiO is a p-type, metal- 
deficient, semiconductor and that the transport and growth properties of NiO 
on nickel are controlled by cationic diffusion of nickel through the NiO lat- 
tice in the presence of nickel ion vacancies created at the gas-oxide inter- 
face.^"^ At present, however, there is considerable controversy as to 
whether the cation vacancies in NiO are predominantly singly ionized^>*^> ^ 
or doubly ionized^ under ordinary conditions of pressure and in the 
temperature range of about 800 to 1700"K. Furthermore, experimental values 
of activation energies for the various kinetic and transport processes are 
not all in accord. 

No systematic study has been performed on the effect of dissociated 
oxygen on the magnitude of the oxidation rate of nickel and on the activation 
energy for this process. Nickel is the base metal in many high-temperature 
alloys for which dissociated gas exposure may occur, in service. We have 
therefore measured the rates of reaction of nickel with atpmic oxygen and, by 
comparing the results with those for the molecular oxygen reaction, have 
developed a model for the mechanism of the process. 


EXPERIMENTAL 


The sample was suspended from an Ainsworth Type 14 recording semimicro- 
balance (10 yg sensitivity) . in a vacuum system and the amount of oxidation 
per unit time was determined from the weight increase. A small bifilar wound 
furnace was placed around the sample inside the vacuum system and was, in 
turn, surrounded by a S-cm^diameter quartz tube connected to the metal vacuum 
system by means of a quick disconnect 0-ring fitting. Figure 1 is a photo- 
graph of the furnace. The core and outer radiation shield were fabricated' 
from boron nitride. The grooved core was wrapped snugly with 0.38-mm (15 mil) 
diameter Pt-13% Rh wire at eight turns per centimeter . (four, turns per centi- 
meter in each direction) and then coated with a ceramic adhesive to improve 
the thermal contact between the wire and the core. The total room temperature 
resistance of the windings was 6.7 fi. The steady-state temperature (in vacuum) 
is plotted as a function of applied power in figure 2; At an inside furnace 
temperature of 1475°K, the winding temperature was 1550°K. The middle 1.5-cm 
section of the furnace was uniform to about 5° at 1175°K. (The sample length 
was 1.2-cm.) 
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The nickel samples were mechanically cleaned by abrasion with 500-grit 
silicon carbide polishing paper, wiped clean,, and washed two cycles with water, 
carbon tetrachloride, and isopropyl alcohol. A specimen was vacuum annealed 
at the beginning temperature of a series of measurements . An initial oxida- 
tion was allowed to occur until the oxide skin exceeded 1000 nm (see Results). 
This ensured also that residual organic surface contaminants were no longer 
present. Subsequent oxidation rate measurements were made on a single sample 
at several different temperatures, usually with cooling between runs. As 
detailed in table I, oxidation experiments were conducted on nickel specimens 


TABLE I.- IMPURITY ANALYSIS OF NICKEL SPECIMENS 


Element 

Concentration 

(ppm) 

Element 

Concentration 

Material 1* Material 2** 

Material 1* 

Material 2** 

Ag 

0.25 


N 


3 

A1 

31 


Na 

0.26 


As 

5.9 


Nb 

0.53 


B 

15 


0 


18 

C 

(Not reported) 

37 

P 

<1 


Ca 

1.8 


Pb 

9.1 


Cd 

<1 


S 

45 


Cl 

2.7 


Sb 

5.9 


Co 

1500' 


Si 

650 


Cr 

290 

2 

Sn 

1.1 

3 . ' 

Cu 

3400 


Ta 

1.2 


Fe 

2500 

12 

Ti 

. 360 


Ga 

. 1.7 


Th 

■ <1 . 


Ge 

0.4 


V 

7 


Hf 

2 


W 

10 ' 


K 

1.3 


Zn 

1.1 


Mg 

400 


Zr 

0.82 


Mn 

6500 


Ni 

98.4% , 

99.99+% 

Mo 

37 






*A11 elements not reported, <0.1 ppm.. 
***A11 elements not reported, <1 ppm. 


of two different purities. Spark source mass spectroscopic, analysis was used 
to determine the impurity level in the less pure material while the vendor 
(Alfa Inorganics) supplied the analysis on the other. The effect of the 
various impurities on the results is considered in the Discussion. 

Oxygen was flowed through the apparatus at a linear flow velocity in the 
range of 40 to 1500 cm sec"^. A microwave cavity of foreshortened 1/4 wave 
coaxial electrical configuration with resonance frequency adjustment (similar 
to cavity 5 of Fehsenfeld et al.^^ but modified to accommodate a 5-cm-diametef 
tube) was placed around the quartz tube and positioned at the center of the 
furnace and sample. It was attached via a waveguide to a Raytheon PGM 2450 
MHz microwave power supply (100 W maximum) which had been modified^^ to improve 
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its stability to 0.1%. With the microwave power supply on, a discharge was 
initiated in the flowing O 2 with a Tesla coil. A stable discharge, which 
partially dissociated the O 2 , could be maintained in the total pressure range 
of 1 to 50 Nm"^ (0.01 to 0.4 Torr) . 

The discharge, as well as the flowing gas, affected the balance readings 
slightly but, as we were interested only in changes of weight, the constant 
extraneous deflections cancelled out. 

The furnace temperature was allowed to reach steady state with the sam- 
ple suspended in vacuo and then O 2 flow was started and the discharge begun. 

The temperature was measured by a calibrated Pt-10%Rh thermocouple 
placed near the sample. The thermocouple junction was covered with a boron 
nitride cap to protect it from the effects of recombining atomic oxygen. Its 
temperature readings were calibrated to those of the sample by comparing its 
readings with those of a 0.25-mm (10-mil) thermocouple imbedded in the sample 
in a series of temperature calibration experiments. The true temperature of 
the sample was obtained by turning off the microwave power supply and extrap- 
olating the recorded temperature back to the instant of power turnoff. In 
this way, extraneous voltages resulting from the microwave field and discharge 
were eliminated. In some later experiments, 25-ym (1-mil) lead wires were 
attached to a 0.25-mm thermocouple imbedded in the sample and annealed. The 
sample temperature was measured, as described above, directly following the 
weight change measurements.. 

/ 

The partial pressure of 0 was determined by titrating with NO 2 at the 
end of an experiment according to the method of Kaufman. ^ ^ In this method, 
the fast reaction / 

I 

NO 2 + 0 NO + O 2 
is followed by the light-emitting slow reaction: 

NO + 0 -»■ NO 2 + hv , 

which is used as an indicator for the first reaction. Accordingly, the light 
output reaches a maximum when the NO 2 flow rate is half that of the 0 flow 
rate and then falls symmetrically to zero at the endpoint at which the NO 2 
flow equals the 0 flow. 

The O 2 flow rate was determined by measuring, near atmospheric pressure, 
the pressure drop across a capillary with a differential manometer^ ^ filled 
with H 2 O. The manometer ha,d been calibrated by observing the flow rate 
through a burette using soap bubbles at a known temperature and pressure. Gas 
pressure at the supply tank was regulated to just aboye ambient pressure by a 
sensitive, 5-psig, regulator (Matheson Model' No. 70) and, following the flow- 
measuring apparatus, was reduced to the system pressure by means of a very 
fine needle valve. 



Nitrogen dioxide flow rates were determined by observing the pressure, as 
a function of time, in a storage vessel of known volume (about 2 liters at an 
initial pressure of approximately 1.5x10** Nm"^) and the temperature as the gas 
was introduced to the system through a needle valve. Flow rates were corrected 
for the amount of association to N 2 O 4 calculated for the average pressure of a 
given flow rate measurement. 

Light intensity measurements were made downstream from NO 2 introduction 
with- an RCA 1P21 photomultiplier equipped with a Corning CS3-72 light filter. 
This combination gave the best filtering of ambient light and good response to 
the emitted light. In practice, we found that the endpoint could be conven- 
iently detected visually in a darkened room just as accurately as by plotting 
the photomultiplier output vs. flow rate. 

The fraction of O 2 dissociation determined by this method was multiplied 
by twice the O 2 pressure, measured prior to turning on the discharge, to obtain 
the partial pressure of 0 . 

We obtained very large fractions of dissociation of O 2 , up to 50%. The 
flow rate was found to be the most important parameter in the control of per- 
cent dissociation of the O 2 . The effects of flow rate, pressure, microwave 
power, and temperature on the percent dissociation of O 2 are shown in fig- 
ures 3, 4, 5, and 6 . The introduction of a metal sample had only a slight 
effect on the measured concentration of atomic oxygen. We could control the 
presence or absence of a "tail" in the discharge in the flowing gases below 
the furnace and its presence increased the amount of dissociation by a fa.ctor 
of 1 . 6 . 

The point of introduction of the N 02 > whether made adjacent to the sample 
inside the furnace or above (downstream from) the discharge zone, had no 
effect on the measured atom concentrations at room temperature and 10 Nm“^ 0 
(35% dissociation) . However, at 1160®K, the measured quantity outside the 
discharge zone was 25% smaller than inside the discharge zone under conditions 
of Pq ~ 1 Nm”2 (0.01 Torr) at 4% dissociation. 


RESULTS 


Past investigators^^” 2 /* of nickel oxidation have shown that a growth rate 
transition occurs at oxide thicknesses of 100 to 1000 nm. At, greater thick- 
nesses, the rate of growth of the oxide is generally found to, be inversely 
proportional to its thickness, x, that is, 

VI 

dx p 

dt X 

In integrated form, then, 

x2 = x2 + 2k't 

op . . 
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or 


+ k t 

op 

where w - w is the total weight increase from t =0. This parabolic rate 
law has been interpreted in terms of the concentration of nickel ion vacancies 
formed at the gas-oxide interface controlling the diffusion rate of cations 
from the metal-oxide interface to the gas-oxide surface. Our measure- 
ments in atomic oxygen were all done with oxide film thicknesses greater than 
1000 nm, and the results were found to follow the parabolic rate law. Values 
of the parabolic rate constants for atomic oxygen were calculated from the 
total rate data by applying small corrections (<5%) for molecular oxygen con- 
tributions. An Arrhenius plot of the atomic oxygen parabolic rate constants 
for nickel is shown in figure 7 for the two materials examined, the data were 
obtained in partial pressures of atomic oxygen varying’ from 1 to 40 Nm"^ (0.01 
to 0.3 Torr) . The parabolic rate constants were independent of atomic oxygen 
partial pressure in this range (see Discussion) . For comparison, the data of 
Douglass et al.^® for air reaction at ambient pressure with nickel, which may 
be taken as representative of molecular oxygen results, were extrapolated to 
7 Nm“^ on the basis of a dependence (as discussed later) and are also 

shown in the figure. The only prior work with the Ni-atomic oxygen system 
was by Dickens et al.,^^ who measured the rates at two much lower temperatures 
and at much smaller reported partial pressures (1 mN m"^ compared to the 
present 10 Nm"^).^ Their data are included' in figure 7; the linear least- 
squares fit to our data for material, for which the most complete and precise 
results were obtained (S.n kp = -11.38 - 13410/T; kp in; g^cm’^sec"^) extrapol- 
ated to their temperatures shows good agreement between the two investigations 

A value of 27 kcal: raole"^ is obtained from the slope of log kp vs. T"^ 
for atomic oxygen reaction with nickel.^ Based on the standard error of esti- 
mate, we estimate the error range to be ±3 kcal mole"-^. 


^hese authors estimated their partial pressures on the assumption that 
enhanced atom , recombination at the higher temperatures of their experiments 
decreased the atom concentration by as much as a factor of 600 from the value 
they measured at room temperature. However, in our apparatus, we observed 
little effect of temperature on the percent dissociation of ©£. It is pos- 
sible, therefore, that the partial pressures of Dickens et al. may be as 
high as 100 times their reported values. This is, however, difficul.t to 
assess accurately because we continuously created the atoms in the furnace 
zone while Dickens dissociated the oxygen upstream of their furnace and 
sample. |: 

^To facilitate comparisons of the present work' with numerous previous 
studies, calories are used as the unit of heat and energy (I'calorie = 4.184 
Joules) . 
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DISCUSSION 

Role o£ Vacancies in Nickel Oxidation 


The parabolic corrosion rate o£ nickel in molecular oxygen is considered 
to be controlled by the di££usion rate o£ nickel ions through the NiO lattice, 
which, in turn, is proportional to the concentration gradient o£ the cation 
vacancies between the gas-oxide inter£ace and the metal-oxide inter£ace. 
Vacancies are created by a series o£ reactions which may be represented by the 
net equations ; 

Ni2++ 1/2 O 2 = 02- + Ni3+ + NiD+ (a) 


NiD+ + Ni2+ = Ni3+ + NiD 


(b) 


2Ni2+ + 1/2 O 2 = 2Ni3+ + NiD+ O^' 


(c) 


The charged species represent ions in the NiO lattice. NiQ^ and NiQ repre- 
sent singly ionized and doubly ionized vacancies for Ni^'*’. It has been 
established^®* that non-ionized vacancies, NiQ , are not significant as 
defects. Because of the small dissociation pressure of NiO, the vacancy defect 
concentration formed at the metal-oxide interface is practically zero and thus 
the vacancy concentration gradient in the oxide lattice between. the metal- 
oxide and oxide-gas interfaces can be taken as equal to the vacancy concentra- 
tion at the gas-oxide surface. 


By a variety of methods, investigators have attempted to determine 
whether reaction (a) or reaction (c) predominates under the conditions of 
their experiments. By applying the mass action law and noting that, for 
reaction (a), [Ni®"*"] = [NiD"’"] (when reaction (b) occurs insignificantly), or, 
for reaction (c) , [Ni®"^] = [2Ni[I]] (when doubly ionized vacancies predominate), 
one can show that, for conditions where an equilibrium concentration of vacancy 
defects exists, for predominately singly ionized vacancies. 


or. 


[NiD^] = K 

■ ; =:C 


1 / 21/4 _ 
a O 2 


P_ exp 
a O 2 ^ 



for the case of doubly ionized vacancies predominating. 


( 1 ) 


[NiQ] 


,:l/3,l/3pU6 , ,-l/3pV6 
C O2 O2 




( 2 ) 


Mole fractions are used for concentrations and [Ni^'*’] = [0^“] = 1 when the 
defect concentrations are small. Note that reaction (a) should produ^e^a 
PO 2 dependence for the equilibrium vacancy concentrations, vyhile a P 02 
dependence is expected for predominately doubly ionized vacancies (reaction (c)) 
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Many investigations (e.g., refs. 1, 2, 4, 6-10), both experimental and 
theoretical, have been performed to detect or predict the relative amounts of 
the two types of cation vacancies in NiO — unfortunately, conclusions from 
these works can support the predominance of either type. In balance, however, 
one factor gives more credence to those proponents of the doubly ionized 
vacancies. Aliovalent impurities in nickel are predominately of valence 
greater than 2, which would tend to increase the pressure dependence of 
vacancy-related properties . ^2 ^nd since, even in high purity nickel, the 
impurity concentration is of the same magnitude as the defect concentration, 
those higher pressure dependence results, i.e., could arise from an impur- 

ity effect. 


Reaction Enthalpies in Nickel Oxidation 

In the following discussion, with the aid of the results from the present 
investigation in atomic oxygen, a detailed examination will be made of the 
various types of cation vacancies which can be present in NiO and control the 
oxidation behavior. 

If the rate of reaction of Ni with O 2 is proportional to the diffusion 
rate of nickel ions in the oxide layer, which in turn is proportional to the 
cationic vacancy concentration formed at the gas/oxide interface, then 

Rate O 2 = C"[NiD] • = :C"[NiD]D^ exp(-AH* . ^^/RT) 

for doubly ionized vacancies. If the vacancy defect concentration is the 
equilibrium value at all temperatures of measurement, then substituting equa- 
tion (2) for [NiQ]: 


Rate O 2 


C"C D P ' exp 
c o O 2 


(l/3AH° (d) + AH* . 


RT 




□Ni is the diffusion coefficient of nickel ion vacancies through NiO (or of 
Ni vacancies, an equivalent process); C", Dq and are constants; AH^iff is 
the net activation enthalpy for Ni diffusing via vacancies in NiO; and 
AHf(d) = AH° is the thermodynamic standard enthalpy of, formation of doubly 
ionized vacancies from O 2 , i.e., the standard enthalpy change for reaction 
(c) . (d and s are used following AH^ to indicate doubly and singly ionized 
vacancies, respectively.) 

Furthermore, 


Rate O 2 = c' 


exp 



RT 



where AH_ is 
Thus, 
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the measured activation enthalpy for reaction of Ni with O 2 



(3) 





(d) 


Similarly, if one assumes that singly ionized vacancies are instead the 
important defects, the following expression can be derived: 

“S, ' “diff * 

As developed below, AH^j^££ may be calculated from experimental values for the 
other terms in equations (3) and (4) . 

^Hf(y) has been determined by several means as tabulated in table II, As 
shown, AH£(y) calculated from diffusion data, using the equation given by 
Mitoff,® is considerably lower (and, in one case, negative) than the other 
values . 


TABLE II.- ENTHALPY OF VACANCY FORMATION IN NiO 


AH^(y), Real mole"® 

Method 

Reference 

19.8 

Gravimetric 

32 


24 ±8 

Gravimetric 

8 


17.8 ±3 

Electrical conductivity 

8 


11.4* 

Diffusion 

8 from 33, 

34 

5.4* 

Diffusion 

calc, from 

35 

-6.4* 

Diffusion 

calc, 'from 

36 

19.1 ±0.7 

Electrical conductivity 

10 


19 

Electrical conductivity 

• Z1 


19.5 ±2 

Weighted mean (neglect- 




ing * values) 




The mean (19.5 ±2 Real mole"^, which neglects the diffusion results and 
is weighted according to the precision of the others) is apparently well rep- 
resentative of AHf(y) . This is also consistent with the observations of Koel 
and Ceilings®® who report a value of 20 Real mole"^ for 1/2 AH£(s) and 
21.7 Real mole"^ for 1/3 AH£(d). (If the vacancies are doubly ionized, then 
1/3 AH£(d) = AH£(y) and,, if they are singly ionized, 1/2 AH^(s) = AH°(y).) 

Reported values for the activation enthalpy for nicRel oxidation by 
molecular oxygen, AHq^, ranges‘*» ® ® from 38 to 68 Real mole”'^. An 

average of all results is 46 ±6 Real mole"®. 


Real mole"® into 


By substituting the values AH 02 = 46 ±6 Real mole"® and 
1/3 AHf(d) = (19.5 ±2) Real mole"® or 1/2 AH£(s) = (19-5 ±2), 
equation (3) or (4), respectively, one obtains a calculated .yalue^ of ^H^j£££ 

of (26.5 ±7) Real mole"®. This value is used in later calculations. 

‘-Recently, Koel and Gellings30 reported an experimentally determined value 
of 24.6 Real mole"® for this quantity. In view of the close' agreement with 
the calculated value, the error ranges given for ^H^j£££ and AH* are' probably 
overestimated. ' ^ 
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Development of Vacancy Models 

We now consider the theoretical relationships between O 2 activation enthal- 
pies and 0 activation enthalpies for. reaction with nickel for the various pos- 
sible combinations of equilibrium vacancy situations where only one type of 
vacancy exists in each environment. These are (1) singly ionized vacancies in 
both cases, (2) doubly ionized vacancies in both cases, (3) doubly ionized 
vacancies in O 2 and singly ionized vacancies in 0, (4) singly ionized vacancies 
in O 2 and doubly ionized vacancies in 0,' (5) doubly ' ionized vacancies in O 2 
and non-ionized vacancies in 0, and (6) singly ionized vacancies in O 2 and non- 
ionized vacancies in 0. Non-ionized vacancies in O 2 have been demonstrated^® 
not to be important. This is also demonstrated by the fact that electrical 
conductivity measurements have shown NiO in O 2 to be a p-type conductor in 
which positive holes (Ni®'*’) are formed which requires formation of ionized 
vacancies.® A seventh situation is possible in 0 — the vacancy concentrations 
(for any degree of ionization) are constant at the same value for all pressures 
and temperatures investigated. 

Equations are tabulated below with the respective AH° symbols and rela- 
tionships following each for the six different equilibrium cases. O 2 or 0 
have been added to the symbols for AH§(d) and AH^(s) to distinguish between 
formation from O 2 or 0, respectively; n following AH° denotes non-ionized 
vacancies . : ■ ' ; 

Case (1) (O 2 single, 0 single): ■ 

Ni2+ + y O 2 = Q2- + + NiO"^ AH°(s,0^) 

0 = y,02 .. AH° = -AH°(0) 

Ni^'" + 0 = Q2- + Ni®-^ + NiD"^ AH°(s,0) = AH°Cs,P 2 ) - AH°(0) 

Case (2) (O 2 double, 0 double): 

2Ni2+ + |- O 2 = Q2- + 2Ni3-^ + NiQ AH°(d,02) 

0 = I O 2 .. . - ; -AH°(0) 

2Ni2"^ + 0 = Q2- + 2Ni3'^ + NiQ AH°(d,0) = AH°(d,02) - AH°(0) 

Case (3) (O 2 double, 0 single): 

2Ni2+ + -| O 2 = 02- + 2Ni®'^ + NiD AH°(d,02) 

Ni®"^ + NiD = Ni^"^ + NiD'^ -AH°(s-^d) 

0 = i O 2 -AH°(0) 

Ni^'^ + 0 = Q2- + Ni®'^ +.NiD'^ AHf(s,0) = AHf(d,02) - AM°(s ^ d) -AHf(O) 
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• 



- ■ 

Case (4) (O 2 single, 0 double) : 


- 


+ 1 O 2 = Q2- + NiS'^ + Nin"^ 

AH°Cs,02) 



Ni^'*' + NiD^ = Ni^"*" + NiQ 

AH°Cs ^ d) 



0 = |02 

-AH°(0) 



'2Ni^'*’ + 0 = Q2- + 2Ni3+ + NiQ AH°(d,0) = AH°(s,02) + AH°(s^d) - AH°(0) 
Case (5) (O 2 double, 0 non-ionized) : 

2Ni2+ + i O 2 = q 2- + 2Ni3‘^ + NiD AH°Cd,02) 

+ NiD= Ni2+ + NiQ"^ .-AH°Cs -> d) 

Ni^'*' + NiO^ = + NiQ^^ -AH°(n ->■ s) 

0 = j O 2 V : 

0 = Q2- + NiD'^'^ AH°(n,0) =,AH°(d,02) - AH°(s->d) 

-AH°(n^ s) - AH°(0) 

Case (6) (O 2 single, 0 non-ionized): 

Ni^-^ + i O 2 = q 2- + Ni3'*''+ NiD'^ AH°(s,G2) 

Ni3+ + NiD"^ = Ni2+ + NiD‘^'^ , -AH°(n -> s) 

0 = j O 2 -AH°(0) 

0 = Q2- + NiD'^'^ AH°(n,0) = AH°(s,02) - AH° (n^- s) - AH°(0) 

Now, by methods similar to those used in deriving equation (3), one can show 
for reaction with atomic oxygen, considering the formation of' doubly, singly, 
and non-ionized vacancies, respectively. 


^^0 ^“diff 3 


(5) 

= ^H*.^f - 1 ^H°(s,0) 


(6) 

AH* = AH*.^^ + AH°(n,0) 


(7) 


For the nonequilibrium situation of fixed vacancy concentration case (7), 
when the temperature is changed, the observed variation in oxidation rate 
would not reflect a change in vacancy concentrations but would be due only to 
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a change in diffusion rate resulting from the temperature effect of the Ni 
diffusion coefficient. Then, 


^“diff 


( 8 ) 


By substituting the appropriate expressions above for H^Cd,0), H£(s,0), and 
Hf(n,0) into equations (5) , (6), and (7), one then obtains the final expres- 
sions for the calculated activation enthalpies for atomic oxygen reaction with 
Ni for the various possible equilibrium situations: 


Case 

(1): 

•k 

* 



AHq = 

^“diff 

Case 

(2): 

* 

k 



AHq = 

^^diff 

Case 

(3): 

k k 

^^0 " ^“diff 

- 1 AH°( 

Case 

(4): 


-1 AH°( 

Case 

(5): 



■ 


% = t AH°(d,02) 

Case 

(6): 

1 

* ★ 

n 


1 ...o, 


1 ...o. 


1 ...O, 


1 ...O. 


1 ...O, 


1 ...b, 


if?' 

\ 1 

1 ...O, 


AHq = + AHJ(s, 02) - AH (n s) - AH^(0| 

And, by equation (8), for the nonequilibrium situation: 

Case (7): 


% = "<iff 


Examination of Vacancy Models 


To determine which of the possible cases are consistent with the experi- 
mental data and which may be ruled out, the following enthalpies are used: 

* , 

AHdiff = (26.5 ±8) kcal mole”^ 

AH^(d,02) = (58.5 ±6) kcal mole"^' 

AH^(s, 02) = (39 ±4) kcal mole"^ ' . 
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= 60.6 kcal mole"^ 

The first three enthalpies were calculated above, while AHf (0) is from the 
JANAF Thermochemical Tables. AH°(s ^ d) = 25 kcal mble"^ (Koel and 
Ceilings^®). Bauer et al.^® estimated that AH°(n t s) = 1/4 AH°(s d) , so 
AH°(n s) = 6.1 kcal rnole"^. 

Our observed value for AHq is 27 kcal mole“^. The calculated values for 
the various cases are summarized in table III. 

TABLE III.- CALCULATED ACTIVATION ENTHALPIES FOR NiO REACTION 


* * 1 
Case Vacancy types AH» (calc.), kcal mole“^ 

O 2 0 


1 

s 

s 

15.7 

±6 

2 

d 

d 

25.8 

±6 

3 

d 

s 

13.0 

±7 

4 

s 

d 

27.6 

±5 

5 

d 

n 

-6.5 

±11 

6 

s 

n 

. -1.2 

± 8.5 : 

7 

- 

fixed 

26.5 

±8 


*d, s, and n refer to doubly-ionized, singlyrionized, and non- 
ionized vacancies, respectively. 


The error analysis of equilibrium cases 1 and 2 is made by noting that 
the first two terms on the right-hand side of the enthalpy equation is, for 
each case, AHq, = (46 ±6) kcal mole“^. For the other equilibrium reactions, 
the error, ±6 kcal mole“S is multiplied by the ratio of the sums of the values 
of the first two terms on the right-hand side of each enthalpy equation to 46. 
While this may not be strictly correct, it is considered a better estimate of 
error range than strictly adding individual uncertainties since the values are 
not independent. 

By comparison of the calculated values of AHq (table III) with the observed 
value (27 kcal mole"*^), we can immediately rule out the existence of cases 1, 

3, 5, and 6. One cannot distinguish, by this analysis, between cases 2, 4, and 
7, which are all nearly equal in value. 


Saturated Defect Model 

Strong support for case 7 is provided by three experimental results from 
the present study, discussed below: (1) the magnitude of the 0 reaction rates 

relative to O 2 rates, (2) the pressure dependence of the reaction rate, and 
(3) the effect of impurities on the reaction rate. 

Comparative reaction fates in 0 and O 2 — For atomic oxygen, one may derive 
(in a fashion similar to that, used for eqs. (1) and (2) for vacancy concentra- 
tions in O 2 ) , for the formation of doubly ionized vacancies, 

[NiDlp '= 4 "^^^ exp[-AG°(d,0)/3RT] (8) 
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Combining equation (8) with (2) and, noting that AG° = AGf(d,02), and that the 
rate is equal to the same constant times the vacancy concentration for O 2 or 
for 0, one may write, for case 2, 


Rate 0 
Rate O 2 


[NiD]Q 


pl/3 

exp[AG°(0)/3RT] 

% 


( 9 ) 


where AGfCO), the standard free energy of formation of atomic oxygen, is 
AG£(d,02) - AG£(d,0). This expression was originally derived by Dickens 
et al.^^ and allows one to calculate a theoretical ratio of rates. 


For case 4, we derive the expression: 


Rate 0 
Rate O 2 


[NiD] 

[NiD'"] 


0 


O2 


4-l/3pl/3 

0 

pi/** 

O 2 


exp 


i AG^Cs.O^), -^AG°(d,0) 
RT 


( 10 ) 


The exponential term in this case can only be estimated. Gulbranson and 
Andrew^^ have calculated a value for ASf(d,02) of -7.5 cal deg"^ mole"^. By 
adding AS° for^the process 0(g) 1/2 02(gy, -15.66 cal deg"^ mole"^, we obtain 

a value for ASf(d,0) of -23.2 cal deg"^ rnole"^. For the formation of a singly 
ionized cation vacancy^ the calculations of Gulbranson and Andrew^ ^ may be used 
to yield a value of ASf(s,02) = -7.3 cal deg"^ mole"^. From the earlier 
derived values for AHf(s,02) (39 ±4 kcal mole"^) and for AH£(d,0) (-2.1 ±6 
kcal rnole"^), we may then write 


AG°(s, 02) = (39,000 ±4,000) + 7.3 T, cal mole"^ 


( 11 ) 


and 

AG°(d,0) = (-2,100 ±6,000) + 23.2 T, cal mole"! (12) 

Note that when equations (11) and (12) are combined for use in equa- 
tion (10), the errors in AHf(s,02) and AH£(d,0) are not additive and, in fact, 
nearly subtract out since they were derived from the same source. A maximum 
combined error is estimated to be ±1 kcal mole"!. The entropy errors, like- 
wise, nearly cancel. The estimated individual entropy errors are ±0.5 cal 
deg"! mole"! and the combined errors less than ±0.1 cal deg"! mole"!, phis is 
equivalent to 0.1 kcal mole"! error in G at llOO^K and so is neglected rela- 
tive to the AH° errors in the calculations . 

A comparison of the theoretical ratio of oxidation rates in atomic and 
and molecular oxygen with that observed is given in table IV for the two pos- 
sible equilibrium cases, 2 and 4. These ratios were determined for the tem- 
perature extremes used in the present experiment. For observed fates in 
molecular oxygen, the data of Douglass et al.^® were averaged with the high 
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TABLE IV.- THEORETICAL AND OBSERVED RATIO OF RATES 




Rate 0 
Rate O 2 


Theoretical ratio 

T,°K 

Theoretical 
Case 2 Case 4* 

Observed 

Observed ratio 
Case 2 • Case 4* 

1050 

67 


30 

2.24 4.30 (2.7 to 7.0) 

1250 

14.4 


6.25 

2.30 4.42 (3.0 to 6.6) 


*The ranges are those resulting from the estimated ±1 kcal mole 
uncertainty in the combined free energies of equation (10) . 


purity data reviewed by Phillips The value of 46 kcal mole"^ determined 
for AHq 2 by Douglass et al.^® is the mean of the values reported in 38-44. 

The atomic oxygen data for both materials represented in figure 7 were aver- 
aged to give representative observed values for that environment. 

From this analysis, one can say, with fair certainty, even in light of 
the estimated theoretical rates, that case 4 (i.e., singly ionized vacancies 
in O 2 and doubly ionized vacancies in 0, both at equilibrium concentrations), 
which differs significantly from the observed rates, does not appear plausible. 
The results of the analysis for case 2 (doubly ionized vacancies for both 
environments) ,' showing a factor of 2.3 variance from the observed, might indi- 
cate that this — case 2 — mechanism is not very probable either. However, in 
view of the range**^* of measured oxidation rates of pure Ni in O 2 by different 
investigators and the observed importances of sample preparation (surface 
preparation, annealing^ pre-oxidation^^ » on the rates, one 
must admit that this lack of agreement between observed and theoretical ratios 
does not completely rule out case 2, and, in fact, case 4 may be marginal. 

The evidence considered thus far, viz, theoretical activation enthalpies 
and rates in atomic oxygen, compared with observed values, within the qualifi- 
cations given in the preceding paragraph, indicates that, while the oxidation 
of nickel apparently occurs under the influence of near equilibrium concentra- 
tions of cation vacancies, in O 2 , the controlling concentration of vacancies 
in the presence of atomic oxygen is less (by a factor of 2 or more than the 
expected equilibrium concentration. Accepting this observation, we hypothesize 
that the atomic oxygen, at the pressure and temperature ; conditions of our^ 
experiments, with its high thermodynamic potential, causes a saturation of the 
surface with vacancy defects (i.e., all defect sites are filled) before an 
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equilibrium concentration can be attained. This model is consistent with the 
suggestion that the vacancy concentration is constant in our experiments 
(case 7 ).d 

Pressure dependence— In support of the proposed saturated defect model, 
we observed no dependence of the oxidation rate of nickel on the partial pres- 
sure of atomic oxygen under the following inclusive test conditions. First, 
all the material 1 specimen points in figure 7, which were obtained at partial 
pressures of atomic oxygen varying from 10 to 16 Nm"^, were extrapolated to a 
common temperature using the experimentally determined activation enthalpy and 
the resultant kp values plotted as a function of Pq (see fig. 8) . There was 
completely random variation of the points in this plot about an average hori- 
zontal line which indicated an average error of ±15% and a maximum random error 
of 29%. These data were obtained for the Ni sample of low purity. Similarly, 
a study was performed on the high purity Ni, material 2, in which the partial 
pressure of atomic oxygen varied between 1 and 45 Nm"^ (thus, pV^ varied by a 
factor of 3.6) and again the variation of the temperature normalized kp's fell 
within experimental error (in this case ±20% maximum). And, lastly, if one 
includes the results of Dickens et al.,^^. which coincide with our extrapolated 
values (fig. 7), at their reported pressure of 10“^ Nm“^ (or if the correction 
of footnote a is valid, 10“^ Nm"^) , the oxidation rate is apparently indepen- 
dent of Pq over a range of 10*^ (or 160 if P is 10"^ Nm“^) . If case 2 or 4 were 
in effect, a Pq/^ dependence should have been observed (eq. (8)). (For ref- 
erence, a pl/3 relationship is indicated on the plot.) 

Impurity effect in 0 and O 2 — Perhaps the strongest evidence for the satu- 
rated defect model (i.e. , case 7 comes from observations regarding the influ- 
ence of aliovalent impurities on the oxidation rate in atomic oxygen. Just as 
changes in pressure would be expected to shift the concentration of vacancies 
at the gas-oxide interface for an equilibrium case, via Le Chatelier’s prin- 
ciple, an equilibrium concentration of vacancies should be changed by the intro- 
duction of impurities of valence different from 2 (see, e.g., Kofstad^^) . Metal 


*^A comment is required at this point regarding the near constancy with 
temperature of the theoretical ratio/ observed ratio (T/0) for cases 2 and 4 in 
table II. One might interpret this as indicating that case 2 or 4 is the 
correct model and that the magnitude at T/0 (which should then be unity) is 
within the range of experimental errors in measuring absolute rates. That is, 
a constant value of T/0 would not normally be expected if the vacancy concen- 
trations were saturated and fixed in atomic oxygen while they varied with tem- 
perature in O 2 . However, because AH£(d,0) is practically zero (-2.1 ±6) 
kcal mole“^), these two equilibrium models for 0 would also predict a concen- 
tration of vacancies in atomic oxygen which remained nearly constant with tem- 
perature changes. (A calculation of AHf(d,0) from the T/0 ratios for case 2 
at the two temperatures yields a value of 1.2 kcal rnole".^ . Thus, either an 
equilibrium model with doubly ionized vacancies in atomic oxygen (case 2 or 
4) or the saturated defect model is consistent with the observed constancy 
of T/0. : ■ ■ 
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impurities with valences greater than 2 would be expected to increase the 
equilibrium concentration of cationic vacancies. As observed by Phillips, 
in his review of nickel oxidation data, the oxidation rate in O 2 of pure nickel 
was about 10% that of less pure nickel. Horn^^ also found that small amounts 
of impurities increased the oxidation rate in 0.2 significantly. Apparently, 
the impurity metals are oxidized rapidly and, in the ionized state, can diffuse 
to the NiO gas surface to affect the Ni vacancy concentrations. 

Most of the present data were obtained with a Ni sample (material 1) , 
which had been represented as being of high purity. Analysis (table 1) showed 
it to contain significant amounts of Mn, Cu, Fe, and Co. Horn^^ specifically 
observed that small amounts of Mn and Cu in Ni increased its reaction rate 
with O 2 ; and Fe (as well as Mn) have been particularly observed to generally 
increase the rate.*^^’^^ Cobalt appears to have a smaller effect on the 

rate.‘*2>47 

In view of the very large concentrations of foreign metals in material 1, 
one would predict that the oxidation rate of this material in atomic oxygen 
should be perhaps a factor of 10 or more higher than that of a very pure 
sample if equilibrium vacancy concentrations were predominant in controlling 
the diffusion rate of Ni through NiO. A few measurements were made on a high 
purity sample (_< 20 ppm metal impurities, material 2, table I. These para- 
bolic rate constants are shown in figure 7 from which it is seen that the pure 
sample actually exhibits a slightly higher rate of oxidation by dissociated 
oxygen (approximately 2.5 times faster than the 98.4% sample). This provides 
additional evidence that an equilibrium mechanism is inoperative. If the 
available defect sites of NiO were all occupied in the pure sample, then the 
addition of impurities, which would have the normal tendency to increase their 
concentration, would be ineffective in doing so. The observations regarding 
the impurity effects may thus be taken as strongly supporting a saturated 
defect model. 


Vacancy Calculations for Saturated Defect Model 

In molecular oxygen at ,0.2 atm and at the two temperature limits of the 
present experiments (1050° and 1250°K), the equilibrium concentrations of 
cationic vacancies in NiO is of the order of 10“^ and 5xl0"^ mole fraction,®’ 
respectively. Since the oxidation rate (at a given temperature) is proportional 
to the number of vacancies (as shown above) , the vacancy concentration during 
the reaction of nickel with atomic oxygen may be estimated. Using the mean 
oxidation rates found (fig. 7) for both materials purities in atomic oxygen 
and the mean rate (as calculated earlier) for pure nickel oxidation by molecular 
oxygen, a vacancy concentration, during oxidation under conditions of the 
piesent experiments, of about 5x10"** mole fraction is calculated for both tem- 
peratures. (The constant value was expected for either of the allowed cases, 
case 7 by definition and case 2 and also case 4 by virtue of the fact that 
AH°(d,0) ~ 0.) 

Again, using the defect concentration equations of Mitoff® and of Tripp 
and Tallen,®^ the temperature at which the NiO defect concentration reaches 
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5x10“** in molecular oxygen at Pg^ = 2x10** Nm"^ (0.2 atm - the "normal pres- 
sure for oxidation experiments) is calculated to be about 1775°K. At 10^ Nm"^, 
it would be approximately 1700°K. These are the. temperatures, at the respec- 
tive O 2 partial pressures, above which one would expect the vacancy sites to 
become saturated if, in fact, the saturated defect model (i.e., case 7) is 
correct for the environment of the present studies. Or, conversely, at 1700°K, 
the vacancy concentration and thus the oxidation rate should become constant at 
pressures greater than 10^ Nm“^. 

Bauer et al.^® measured oxidation rates of Ni in O 2 at pressures up to 
2 MNm'^ between 1275°K and 1475°K, At the higher pressures, above about 
5x10^ Nm“^, their observed rates are apparently constant at all temperatures. 
Further, it appears that the pressure at which the rate becomes constant may 
be smaller as the temperature increases. The "break" in their 1475°K data is 
at approximately 5xio? Nm"^. At that temperature and pressure, the vacancy 
concentration is calculated to be 3xl0"**. This is in excellent agreement with 
the value obtained in the present study, particularly in view of the simplify- 
ing assumptions and, more spectacularly, the experimental uncertainties. 

One would expect, however, that the high-pressure results of Bauer et al.^® 
for O 2 should yield an observed activation enthalpy identical to that observed 
in the present study, i.e., it should reflect only the vacancy diffusion acti- 
vation enthalpy. Bauer in fact, did, not observe such a decrease in the acti- 
vation enthalpy in the high-pressure range. This is not consistent with their 
observed pressure independence and their model. Note that their statement 
regarding the closeness of their observed activation enthalpy (50.5 kcal) to 
that for diffusion of Ni in NiO must be modified to recognize that the measured 
activation enthalpy for Ni diffusion also includes the term for vacancy forma- 
tions. Thus, their conclusion regarding the smallness of the enthalpy of 
vacancy formation is incorrect and they should not have ignored the fact that 
their Arrhenius slopes were the same in both the low- and high-pressure regions. 

Realistically, it is possible to consider that their data do not reflect 
a clear independence of pressure in the higher range (fig. 2 of Bauer et al.®®). 
If this is accepted, the minimum saturation vacancy concentration would have 
to be greater than that calculated in O 2 at 1475°K and 2MNm“^, or 3.8x10”**. 

This value is still consistent with that observed in atomic oxygen. 


CONCLUDING REMARKS 


By considering the results of past nickel oxidation (and related) studies 
in molecular oxygen with the present results in atomic oxygen in terms of all 
possible NiO defect structure combinations in the two environments, the fol- 
lowing equilibrium vacancy situations (designated as in Discussion) have been 
definitively ruled out: (1) singly ionized vacancies forming at equilibrium 

in O 2 at ordinary pressures and in 0 in the pressure range of our experiments, 
(3) doubly ionized vacancies in O 2 and singly ionized vacancies in 0, (4) 
vacancies singly ionized in’ O 2 and doubly ionized in 0, (5) doubly ionized 
vacancies in O 2 and neutral (non-ionized vacancies in 0,; and (6) singly ionized 
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vacancies in O 2 with non-ionized vacancies in 0. Except for case 4, these were 
eliminated on the basis of observed activation enthalpies of reaction compared 
with theoretical. The third situation appeared not .too probable since the 
observed oxidation rate in atomic oxygen was much smaller (outside of experi- 
mental and calculated uncertainties) than the theoretical rate. Additionally, 
all the evidence cited below, in favor of the saturated defect case, is unfav- 
orable toward this fourth situation, .as well as the other two equilibrium cases 
above . 

On the basis of present evidence, two situations are possible: (i) doubly 

ionized vacancies formed in the presence of O 2 and also in the environment of 
dissociated oxygen, both at equilibrium concentrations near the surface, and 
(ii) a saturated and constant concentration of vacancies . (of any degree of 
ionization) formed in the presence of atomic oxygen, irrespective of the situ- 
ation for O 2 . These are both definitely possible on the basis of the observed 
activation enthalpy of reaction of Ni with 0. Observed rates in 0 and in O 2 
and, thus, calculated defect concentrations, provide information that can be 
interpreted in terms of either model. But the apparent pressure independence 
of Ni reaction rates with atomic oxygen did indicate that the saturated defect 
model is the correct one. However, a Pq ^ effect, which would have supported 
the equilibrium case, is not large compared to experimental uncertainties and 
therefore the apparent absence of such an effect should not be construed as 
absolute proof for the saturation of defect sites. Finally, the fact that 
large concentrations of Mn, Cu, Fe, and Co present as impurities in the nickel 
did not increase the oxidation rate is quite strong evidence that atomic oxygen 
at the pressures and temperatures studied creates a saturation of nickel in 
vacancy defect sites in the NiO film. 

The activation enthalpy for reaction of atomic oxygen with nickel appears 
to be reasonably established, within the limits of error for the theoretical 
value. The experimental value of 27 kcal mole"^ determined from rate measure- 
ments on the impure sample, and which agree well with the lower temperature 
results of Dickens et al.,^^ may be considered as possibly a minimum value 
when compared with the few data points obtained with the pure material. Yet 
the maximum possible theoretical value is (26.5 ±8) kcal mole"^ and so it is 
reasonable to expect that the true activation enthalpy lies somewhere within 
the range of 27 to 34 kcal mple"^, the lower limit being more probable. 

Additional work should be done, however, on accurate determinations of the 
reaction rate- of pure Ni arid, particularly, its oxygen partial pressure depen- 
dence, in order to more confidently verify the saturated defect model. Specific 
attention must be given to sample preparation methods for comparisons of data. 

Electrical conductivity measurements, which can normally be performed 
more precisely than oxidation rate measurements, should be made to determine 
pressure effects in atomic oxygen. A value for the activation enthalpy for 
conductance, AH^j, in atomic oxygen is predicted to be 5.5 kcal mole"^, the 
same as in molecular oxygen® if the saturated defect model is controlling; or 
AH* (in O 2 ) + 1/3 AHf(d,02) - 1/3 AHf(O) = (4.8 ±2) kcal mole"^ if equilibrium 
doubly ionized vacancies are present in both environments. If singly ionized 
vacancies are equilibrated in O 2 but doubly ionized vacancies in 0 (the only 
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other reasonably possible equilibrium case since others would require a nega- 
tive AH(j in atomic oxygen) , the predicted AH(j for 0 is AH 0 (in O 2 ) + 1/3 
AHf(s,02) + 1/3 AHf(s d) - 1/3 AHf(O) or (6.6 ±2) kcal mole-^. 

Finally, studies in molecular oxygen (oxidation of Ni or conductivity of 
NiO) should be extended to high pressures and temperatures to attempt satura- 
tion of defect sites to test the predictions resulting from the atomic oxygen 
data. Measurements in atomic oxygen should be extended to as low pressures as 
possible so that verified equilibrium concentrations of vacancies can be 
determined. (Temperature should have little effect since AH£(d,0) ; 0.) Then 
it could be shown conclusively which degrees of vacancy ionization occur in O 2 
and in 0. 
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Figure 1.— Disassembled view of reaction furnace. 



Temperature, ®K 




f 





Flow rate O2 , cm^ min~* (STP) 

Figure 3.— Effect of flow rate on percent dissociation of oxygen. Conditions: 

- Ni sample in furnace; 80% microwave power; all discharge in furnace; 
po = 15 Nm-2; T = IO75OK. 
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Effect of microwave power on percent dissociation of oxygen; 100 
power =100 W. Conditions: Ni sample in furnace; all discharge 
in furnace; P£ = 15 Nm~^; flow rate = 8.0 cm^min"^ (STP) ; 

T = 1075°K. 2 
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Figure 6.— Effect of temperature on percent dissociation of oxygen. 

Conditions: 80% microwave power; all discharge in furnace; 

P° = 15 Nm“2; flow rate = 8,0 cm^min"^ (STP) . 
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Figure 7.- Arrhenius plots for nickel reaction with atomic and molecular oxygen 
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Figure 8.— Pressure dependence of parabolic rate constant' for atomic oxygen 
reaction with nickel; all rates extrapolated to 1152°K. 


